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Abstract 



We prove the existence of a weak solution to the three-dimensional steady compressible 
isentropic Navier-Stokes equations in bounded domains for any specific heat ratio 7 > 1. 
Generally speaking, the proof is based on the new weighted estimates of both pressure and 
C^ , kinetic energy for the approximate system which result in some higher integrability of the 

density, and the method of weak convergence. Comparing with [12] where the spatially 
periodic case was studied, here we have to control the additional integral terms of both 
pressure and kinetic energy involving with the points near the boundary which become 
degenerate when the points approach the boundary. Such integral terms are estimated 
using some new techniques, i.e., we use the techniques of the mirror image and boundary 
^^ , straightening to prove that the weighted estimates of both pressure and kinetic energy for 

r^ ' the points near the boundary can be controlled by the weighted estimates for the points 

on the boundary. Moreover, we prove that once the weighted estimates of the kinetic 



C^ ' energy in the direction of the unit normal to the boundary are bounded, we can control 



the weighted estimates of the total energy on the boundary. 

Keywords: Steady compressible Navier-Stokes equations, existence for any 7 > 1, weighted 
estimate, bounded domains, viscous flux. 



1 Introduction 

In this paper we shall prove the existence of a weak solution (p, u) to the following steady 
isentropic compressible Navier-Stokes equations in a bounded domain J7 C M^ for any specific 
heat ratio 7 > 1: 

div(pu) = 0, (1.1) 

-^Au - /iVdivu + div(/9u (S,u)+\7P = pi + g (1.2) 

with Dirichlet boundary condition 

u = on dn. (1.3) 

Here u = {ui,U2,U3) is the velocity and p is the density, the viscosity constants fi and jl 
satisfy fj,>0, fl = p + X with A -|- 2///3 > 0, the pressure P for isentropic flows is given by 

Pip) = ap^ 



with a being a positive constant and 7 > 1 the specific heat ratio, f = (/i,/2,/3) and 
g = (91,02,93) are the external forces, and for simphcity, we assume that 

f, s(^L^m. 

Moreover, the total mass is prescribed: 

pdx = M > 0. 



n 

In the last decades, the well-posedness of the equations (1.1), (1.2) for large f and g has 
been investigated by a number of researchers. In 1998, under the assumption that 7 > 1 in 
two dimensions and 7 > 5/3 in three dimensions. Lions [13] proved the existence of weak 
solutions to boundary value probelms for (1.1), (1.2). Roughly speaking, the condition on 
7 comes from the integrability of the density p in L^(i7). The higher integrability of p has, 
the smaller 7 can be allowed. If f is potential, then weak solutions are shown to exist for 
any 7 > 3/2, see [15]. Then, Frehse, Goj, Steinhauer, Plotnikov and Sokolowsk established 
an improved integrability bound for the density by deriving a new weighted estimate of the 
pressure in [5, 16], where the authors assumed a priori the L^-boundedness of pu^ which, 
unfortunately, was not shown to hold. Later, by combining the L°°-estimate of A^^P with 
the (usual) energy and density bounds, Bfezina and Novotny [1] were able to show the 
existence of spatially periodic weak solutions to (1.1), (1.2) with some symmetries for any 
7 > (3 + \/41)/8 when f is potential, or for any 7 > 1.53 when f G L°°, without assuming 
the boundedness of pu^ in L^. Recently, Frehse, Steinhauer and Weigant [8] established the 
existence of weak solutions to the Dirichlet problem for (1.1), (1.2) in three dimensions for any 
7 > 4/3 in the framework of [1]. Also, the existence of weak solutions in the two-dimensional 
isothermal case (7=1) was obtained in [6, 7]. 

More recently, the authors of this paper [12] proved the existence of a spatially periodic 
weak solution with the same symmetries as in [1] to the system (1.1), (1.2) for any specific 
heat ratio 7 > 1. The proof is based on the weighted estimates of both pressure and kinetic 
energy for the approximate system which result in some higher integrability of the density, 
and the method of weak convergence. 

The aim of this paper is to prove the existence of a weak solution to (1.1)-(1.3) in bounded 
domains in M^ for any specific heat ratio 7 > 1. Comparing with [12], here we have to 
establish the estimates involving with the points near the boundary which are degenerate. 
So, we cannot use the same techniques as in the derivation of the estimates of the interior 
points in [12] to get the weighted estimates of both pressure and kinetic energy for the points 
near the boundary. Instead, such weighted estimates are established using some new test 
functions, see the paragraph below Theorem 1.2 for more details on the proof idea. 

We mention that for a 3D model of steady compressible heat-conducting fiows (i.e., 
the steady compressible Navier-Stokes- Fourier system), Mucha and Pokorny [14] recently 
studied the existence of weak solutions under some assumptions on the pressure and heat- 
conductivity, which unfortunately excludes the case of polytropic idea gases. For the corre- 
sponding non-steady system to (1.1), (1.2) with large initial data. Lions [13] first proved the 
global existence of weak solutions in the case of 7 > 3n/(n + 2) (n = 2, 3: dimension). His 
result has been improved and generalized recently in [2, 3, 4, 10, 11] and among others, where 
the condition 7 > 3/2 is required in three dimensions for general initial data. 

Before defining a weak solution to (1.1)-(1.3), we introduce the notation used throughout 
this paper. 

Notation: Let G be a domain in M^. We denote by U'{G) the Lebesgue spaces with 
norm || • \\lp, by W^'^{G) {k G N) the Sobolev spaces with norm || • ||vk'"'P; by C^{G) (resp. 



C^{G)) the space of /cth-tinies continuously differentiable functions in G (resp. G). In 
particular, F™(G) = VF™'2(G) with || • \\h^. By V'{G) we denote the dual space oiV{G). 
Bji(a) := {x € M^ : |x — a| < R} denotes the open ball centered at a with radius R. The 
capital letter C (sometimes used as C{X, ¥,•••) to emphasize the dependence on X,Y,- • •) 
denotes a generic positive constant which can vary from line to line. 

Now, let us recall the definition of a bounded energy weak solution to (1.1)-(1.3). 

Definition 1.1. (Renormalized bounded energy weak solution) We call (p, u) a renormalized 
bounded energy weak solution to the system (1.1)~(1.3), if the following is satisfied. 

(1) p>0, p£ L^'{n), u G H^{n), J^ p{x)dx = M > 0. 

(2) (/?, u) satisfies the energy inequality: 

/ {p\Vu\^ + fi\dwu\^)dx < {p{ + g)-udx. (1.4) 

Jn Jn 

(3) The system (1.1), (1.2) holds in the sense ofV{Q). 

(4) The mass equation (1.1) holds in the sense of renormalized solutions, i.e., 

dw[b{p)u\ + [b'{p)p - 6(p)]divu = in V'{n) (1.5) 

for any b G C"'^(M), .such that b'{z) = when z is big enough. 

Then, the main result of the current paper reads as 

Theorem 1.2. Let Q. be a G"^ -domain in M^, f, g G L°°{^). Then for any 7 > 1, there exists 
a renormalized bounded energy weak solution (p, u) to the system (1.1)-(1.3). 

Roughly speaking, the proof of Theorem 1.2 is based on the uniform a priori estimates 
for the approximate solutions and the weak convergence method in the framework of Lions 
[13]. Comparing with [12] where the spatially periodic case was studied, here we have to 
control the additional integral terms of both pressure and kinetic energy involving with the 
points near the boundary which become degenerate when the points approach the boundary. 
Such integral terms are estimated using some new techniques, including the techniques of 
the mirror image and boundary straightening. The basic steps of the proof are the following: 
First, we use the distance function as a test function to get the weighted estimates of both 
pressure and kinetic energy in the direction of the unit normal to dVi, i.e. the weighted 
estimates of Pg + Pei'^e ' V(i(x))^ for the points on dO,, where \7d{x) is the normal direction 
of the boundary at xq satisfying d{x) = \x — xo\. Second, by straightening the boundary 
locally, we can show that the weighted estimates of the total kinetic energy for the points on 
do, can be controlled by exploiting the weighted estimates of Ps{ug • V(i(x))^. Then, by using 
the method of the mirror image which is often used to construct the Green's function of the 
first kind for a ball or a half-ball, etc, we can construct a special test function to prove that 
the weighted estimates of both pressure and kinetic energy for the points near the boundary 
can be controlled by the weighted estimates of the points on the boundary dO,, and thus the 
weighted estimates for both pressure and the total kinetic energy can be closed, leading to the 
desired a priori estimates, in particular, the higher integrability of the density. Finally, with 
the help of the established a priori estimates, we can take to the limit in the same manner 
as in [12] to complete the proof of Theorem 1.2. 

This paper is organized as follows. In Section 2, we first construct a sequence of approx- 
imate strong solutions (/9e,U£) and then in Subsection 2.2 we derive the uniform weighted 
estimates for both pressure Pg and kinetic energy yOe|ue|^ in Q.. In Subsection 2.3 we take 
the first level limit e — )■ and Eq ^ and then show the additional uniform estimates for 
the velocity u^ and the pressure Ps in terms of the quantity A = \\P^\us\'^ + P5 [u^l^^^^H^i, 



< a < 1. Then, by using a bootstrap argument, we prove that A is uniformly bounded 
which in turn imphes the uniform //^-boundedness of U5, and the L''-boundedness of Ps, 
psug and p^ju^p for some r > 1. In Section 3, we prove the main theorem by using the weak 
convergence method in the framework of Lions [13]. 

2 Uniform estimates of the approximate solutions 

2.1 The approximate system 

To prove Theorem 1.2, we first work with the standard approximation problem in 0, with 
positive parameters e,d, Eq < 1: 

eo {p-h) + div(pu) - eAp = 0, (2.6) 

eo(/9 + /i)u - ^Au - /iVdivu + div(/9u (g) u) + eVp • Vu + VP^ = pf + g (2.7) 

with boundary conditions 

u = a„/9 = on dVL, (2.8) 

where h = M\Vl\''^ , Ps = p^ + 6p'^, and n is the outer normal vector to d^. 

According to [15], we have the following existence result for the problem (2.6)-(2.8). 

Proposition 2.1. There is at least one strong solution (p£,U£) to the problem (2.6)-(2.8) 
with the following properties (7 = max{7,4}^.- 

Pe(^H'^{^), \ie(^Hl{Vl)r\H'^{Vl), ps>Oa.e.inn, p^dx = M; (2.9) 

Jn 
WPsWh'^ + ||u,||^2 < C{eo,e,6,M,i,g,h); (2.10) 

||u,||^i < C{n,h,{, g){l + llp^llis), (2.11) 

ellVp.lli^ < C(f],/i,f,g)(l + llpellLs)', (2.12) 

S\\p,\\LS<C{n,h,{,g). (2.13) 

Before passing to the limit e — )■ 0, Eq — ^ 0, (5 — )■ to get the existence of a weak solution 
to the system (1.1), (1.2), we need first to show necessary uniform weighted estimates for 
(ps,Us). In fact, we will prove the following theorem: 

Theorem 2.2. Let {pi,, u^) be the solution of the approximate problem (2.6)-(2.8) established 
in Proposition 2.1. Then, for any (5 > 0, there exists a constant 5q > 0, such that for any 
e,eo < ^0, 

I ^^f^-^^dx < C(l + ||P,||i. + \\pe\vie?\\L^ + ||Ue||^0 (2.14) 

JO, F ~ ^0| 

for all a G (0, 1), xq G ^, where the constant C depends on f, g, p, p, M , 7, Vt and a, but 
not on 6, e, Eq and xq. 

The proof of Theorem 2.2 is broken up into several lemmas given in Subsection 2.2. 

2.2 A potential estimate 

In this subsection we will derive both the interior and boundary weighted estimates for P^ 
and /3e|ue|^ which can be understood as estimates in the Morrey space. 



Lemma 2.3. Let (p^jMg) be the solution of (2.6)-(2.8) given in Proposition 2.1. Then the 
following estimate holds. 

I ^^f^^-^^dx < C(l + llP.ll^i + ||p.|u,p||^. + ||u,||2,,) (2.15) 

for all a G (0,1), xq G ^ and Rq = -^dist^xojdQ), where the constant C depends on f, g, 
^, Rq, jl, M , 7 and a, hut not on 5, e and Eq. 

Proof. For xq G il, we define (f) = {4>^ , (fP' ,(f)'^) with 

.^<(.) = l£;3!i:„!£^), ,: = 1,2,3, 

|x — Xol ^0 

where < a < 1, r/ G (7°°(0, oo) is a cut-off function satisfying < ?7(t) < 1, \Drj\ < 2 and 

*) = (' M-' 

[O |t| > 2. 
Testing (2.8) by (p and performing a direct computation similar to that in [12], we obtain 

f , ^'- , dx + (1 - a) /" j^^^^dx 

< c(l + \\Pe\\L^n) + ll/^elUePllLi + ||Ue H^^i + eollPeU^H^i + e^ || Vp,||2 ^ V (2.16) 



If we take eo < (^ and e < 6^, then (2.15) follows from (2.12), (2.13) and (2.16) immediately. 

D 

For (T > 0, let us set T„ = {x & Cl \ d{x) < a}. The following lemma relates the 
smoothness of the distance function d in Tu to that of the boundary dO, (we refer to the 
Appendix of Chapter 14 in [9] for the details). 

Lemma 2.4. Let Q be bounded and 50 G C^ for k > 2. Then there exists a positive constant 
a depending on Q such that d G C (Fq-). 

Now, we are in a position to get the uniform weighted estimates for the points near the 
boundary. The following lemma is inspired by Lemma 3.3 in [8], but the difference lies in 
that we prove in addition a weighted estimate for the kinetic energy in the direction of the 
unit normal to dVt for the points on dVt which is crucial in our proof. 

Lemma 2.5. Let [pi,,Uf,) be the solution of (2.6)-(2.8) given in Proposition 2.1, and o" > 
be a constant such that d{x) G C'^iVa-), Then, the following estimate holds 

P,|Vd(x)P+fe(u,-Vd(x))% ^^,^ , „„„ 2,1 ^||„||2N ..-,7^ 

. dx <C[l + \\Pe\\L^ + \\pe\Ue\ \\l^ + \\Ue\\m) (2.17) 

for all a G (0, 1) and xq G dO., where the constant C depends on f , g, p, /i, M, 7, a and a, 
but not on 5, e, Eq and xq. 

Proof. For xq G dVi, a G (0, 1), m = 2^^, we define 4> = {4>^ , (j)^ , (p^) with 

,,-, , d(x)did(x) /2|x — xol\ ^ 

^\x) = ^,^ / ' ^ ^, ^ ?7 J ^ in [7, i = 1,2,3. 

^ ^ ' (d(x) + |x-xo|™)° '\ a J 



/. 



Multiplying (2.8) with (p and integrating over fi, we find that 

P£div(pdx + / peu\uidj(t)'dx = / (eo(P£ + h)u.ir + eVu^ • Vp£)4>dx 

Jnn J B„{xa)r]n JB^(xo)nn 

+H / V Us : V (pdx + jl / divug div(/)(ix + / {psf + g)(pdx. (2.18) 

J B„(xn)c\Q. JB„(x,))nn JB„(x,))nn 



iBa{xo)nn JB„{xo)nn JB„{xo)nn 

A straightforward computation gives that 

djd{x)did{x) + d{x)dijd{x) ^2\x — xq\ 



d^<t>\x) = ^t:t J ::r ' ^i 



((i(x) + |x-xo|™)" a 

dix)did{x)\djd{x) + dj\x — xqV^] ,, ,, d(x)did(x) „ ,2|x — xqI 

-« .^.^N ,1^ ^ iJ,A»^i ^(1^ - ^ol) + 717-^-1-— -i;;x7: 9M- 



(<i(x) + |x - xo|'")"+i (d(x) + |x-xo|™)" -^ 0- 

d(x)dijd(x) ,2|x — xnL d(x)didix) _ ,2|x — xnL 

Vi 1 ) + 717Z^^—i;;^9jr]{- ^ 



((i(x) + |x-xo|'")" 0- ' ((i(x) + |x-xo|™)" -^ cj 

(1 -a)d(x) + |x-xor„ ,, .r. ,. ^ .2|x -xqU 

+ , ,, / , ■ — ; — - — --rOid[x)djd(x)ri[ ) 

2(d(x) + |x-xo|™)"+^ » V ; J V ;/v ^ j 



(1 — a)d(x) + \x — xqI'^ ,^ ,, , ,„•, ,,,_ ,, , ,^,-, ,, ,2|x — xnl 
+ lT717Z^^T^—^n;W^i9^d{x) - M\x)){d,d{x) - M^(x))7?(- ' °' 



2{d{x) + \x- Xo|™)"+i V * V ; v ;n j v / \ j) i\ ^ 

2a[djd{x)di\x — xo|™ — did{x)dj\x — xqI™] ,2|x — xqI ^ 



2(d(x) + |x-xo|™)"+i cr ^ 

a^d^(x)g»|x - xo\"^dj\x - xpp 2|x - xpL 

^2(d(x) + |x-xoh)°+i((l-a)4a;) + |x-xoh)''^ a ^' ^^ ^ 

where Mi{x) = a(i(x)(9i|x — xo|'"((l — a)(i(x) + |x — xqI™")"^, i = 1,2,3. From (2.19) it follows 
easily that 

div0(x) > (1-«)|V4^)P ,( ^l^-^ol ) _ c. (2.20) 

^^ ^ - 2(d(x) + |x-xo|™)° CJ ^ ^ ^ 

The advantage of the representation (2.19) lies in that one can clearly identify the sym- 
metric and skew-symmetric parts. From (2.19) and (2.20) we get that for any xq G dQ, 



L 



PsA^yHx)dx> f (1 a)P,jVd(x)p _^^^^|| II ^2.21) 



and 



3 



X- f ^AM ^ /" (l-a)pe(u,-Vd(x))^ 

> / PeUiUjdjd) dx> I ^, ,, , ; 1 r 

ij^jBAcco)nn ' '' ' -JB.{xo)nn 2{d{x) + \x - xo\^r 

-C{a)\\p,\u,\^\\Li. (2.22) 

Substituting (2.21) and (2.22) into (2.18), we obtain 

P,|Vd(x)|2 + p,(u,-Vd(x))2 



/. 



-dx 



Ba(xo) |x-xor 

c 

< -, (1 + ll^ellLi + ||p,|u,p||ii + WueWJj^+eoWpsUsWLi+e^WypsWh)- (2.23) 

1 — a ' 

Thus, taking cq < (^ and e < (5^, we see that (2.17) follows from (2.12), (2.13) and (2.23) 
immediately. D 



For any xq £ dQ, since the Navier-Stokes system is invariant under the rotation trans- 
formation, we can assume that the X3 coordinate axis lies in the direction z^(xo), here i^ixo) 
is the unit inner normal to dfl at xq. Consequently, one has V(i(xo) = i^ixo) = (0,0,1). 
Besides, by Lemma 2.4, d £ C^(ro-), ||d||c2(r^) < C and for any e > small enough, whence, 
there exists a constant eo > such that for eg = minjeo, o"}, 

||Vd(x)-(0,0,l)||ci(ij(^„,e„)nr.) < e, 
which combined with Lemma 2.5 implies that for any xq G B{xq, ^) n dQ, 



L 



B(5;o,f)nr<, \x-xq\ 



-dx < C{1 + llP^Jlii + llPelUel'llii + llUell^i) (2.24) 



Now, we define the following transformation 99 which straightens the boundary T = 
B{xq, eo) n 5il in iV = B{xo, eo) n F^ to 

yi=Xl, y2 = X2, y3 = d{xi,X2,X3), 

and set 

yo = ip{xo), N' = tp{N), T' = ip{T), Ue(y) = Ue(x), pe{y) = ps{x). 

Then, under this transformation, the system (2.6)-(2.8) is transformed to (here all derivatives 
are with respect to y unless noticed explicitly): 

eo{pe -h)+ div{peUe) - e^Pe = FliPe, d, Ue), (2.25) 

eoipe + h)us - /uAUj - /iVdivUe + div(peUe Ug) + VPs^ips) + eVps ■ VUe 

= Pe{ + g + F2{pe,d,Ue), (2-26) 

where F2 = {F^,F^,Fl) and 

Fl{ps, d, Ue) = -{dx,d - (53j) • dy^ipeUi^) 

-e[2dx,ddi3Pe + 2dx,dd23pe + d^Pe^xd + (|V^dp - 1)5|/J,], 

F^{pe,d, Ue) = n[2dxiddi3Uj^ + 2da:2dd23Uj, + dsUj^Axd + (iV^dl^ - 1)9|%] 
-il[{dx^d- 53j)d3divue + dx^xidd3Ui^ 
+{dx,d - 63i) ■ {d3jUi^ + {dx^d- S3j)d^Ui^)] 

-{dx,d - 53i) ■ {dysipeUi^Uj^) + e[d3PediUje + diped3U^ + {dx,d - 63i)d3ped3Uj^]} 

-{dxjd- 53j)d3ps, 

= : 0{e)F{V^Ue, Vpe ■ Vu„ d3{peUi,UjJ) + 0{\VUe\). 
If we define e~o := jr^^o, then (2.24) implies that for any yo £ B{yQ, ^) H T' , one has 

/ - ^f^ ^ ^-^ dy < C(l + llP^ll^i + \\pe\ue\^\\L^ + ||u,|||0- (2-27) 

JB{yo,f)nN' \y-yor 

Next, we show that with the help of the estimate (2.27), we can also control the weighted 
estimates of the total kinetic energy for the degenerate points on the boundary. In fact, one 
has the following lemma. 



Lemma 2.6. For any y^ G B{yQ, ^) n T' , we have 

iB{yo,f)nN' \y-yo 



f - —2 

/ _ -^^^dy<C{l+\\Ps^\L^ + \\pe\vie?\W + \\vie\\m). (2.^ 



Proof. Set N' = {(2/1,2/2,2/3)! (2/1,2/2, -2/3) G A^'}, iV' = A^' U N' and 2/' = (2/1,2/2)- Since p,, 
Ue G W^^'^(iV') and u^ = on T, we denote 



Ue(2/',2/3) 



-Ue(2/',-2/3) for 2/3 < 
Ue(2/',2/3) for 2/3 > 0, 



Pe 



Peiy'^-ys) for 2/3 < 

Pe(2/',2/3) for 2/3 > 0. 



to deduce that u^ G H^^'^(iV'), pe G H^-'^'^(iV') and the foHowing equations hold. 

— div(/5eU£ (g) Ug) = Gg, a.e. in N' 

where 



(2.29) 



GK2/',2/3) 



{eoipe + h)us - fiAus - ^Vdivue + VP5, - Pef + eV/>e • Vue - g}(2/', 2/3) 

--^2(Pe,'^,U£)(2/',2/3) for 2/3 > 0, 

{eoipe + h)ue - /uAue - //VdivUe + VP5, - p^f + eVp^ • Vuj - g}(2/', -2/3) 

,--^2(Pe,'^,Ue)(2/',-2/3) + 2«93(Pe%%)(2/', "2/3) for 2/3 < 0. 



For any yo G ^(2/0, f )nr', multiplying (2.29) by ^^(2/) 



Jy-yoY „^ '2\y-yo 
\y-yo\^ 



v{- 



eo 



and integrating 



over N', using the fact that 0*(2/',2/3) = '/'*(2/', ~2/3) = 0*(2/',O) on T', we obtain by a direct 
calculation that 



:i-a) 



+ 



Ps, + PeU 



B(j/o,f) iy-yo 



^d2/ < C(l + llP^JI^i + \\pe\Ue\'\\Li + ||u,|||^i; 



Pe\Ui,U3^ 



s(go,f) 12/ -yo 



^d2/ + eo||/5eU^||Ll+£ l|V/9£||22 



(N')- 



(2.30) 



If we take Eq < 6 and e < 5^, then (2.28) follows from (2.12), (2.13), (2.27), (2.30) and 
Young's inequality immediately. D 

Then, by using the method of the mirror image which is often used to construct the 
Green's function of the first kind for a ball or a half-ball, etc., we can construct a special 
test function to prove that the weighted estimates of both pressure and kinetic energy for 
the degenerate points near the boundary can be controlled by the weighted estimates for the 
points on the boundary: 



Lemma 2.7. For any y G B{yo, ^) n N' , one has 

L 



iB{y,f)nN' \y -y 



Ps^ + p-,u|^^ ^ ^^^ ^ ii^^^ii^^ ^ Wpsln^n^^ + ||u,f^,; 



(2.311 



Proof. For any y = (2/1,2/2,2/3) G -8(2/0, f) n iV', let y* = (2/1,2/2,-2/3), 2/0 = (2/1,2/2,0), and 
define 



i^\y) 



{y-yy (l-2(53.)(2/-r) 



^ .2^/--2/| . 

Vi — z ), « = 1,2,3, 

eo 



\y-y\ \y-yr 

thenipiy) = (V'^V^V'^) = on diB{y,eo) D N'). 

Multiplying (2.26) with ijj and integrating over A^', since \y — 2/0I ^ \y ~ y*\ for all y £ N' 
and 



5,VX2/) 



{y-yy (l - 263,)iy - y*y 



12/ -2/1' 



12/ -2/ 



5j^(l2/ - 2/1) 



+ < 



I2/-27I' 



a 



(2/ - y)*(2/ - yy (1 - 2.53.)<J^, , (1 - 2<53.)(2/ - r )^(2/ - y*)' 



\y-y 



a+2 



\y-y 



+ a- 



\y-y 



*|a+2 



-}^(l2/-2/|). 



we find that 

/ , - , 

lB{y,^)nN' JB{y,^)nN' \y-y[ 

h 



[ . Ps^dwijdy > {3 - a) [ . _^<iy-C(eo)||P5j|Li 

jB(vM)nN' JB(v.^)nN' \y-y\ 



iB{y,^)nN' \y - yo 

and 

3 /• r 7, IV-, |2 



C{h) / _ I.. l.^ dy (2.32) 



r r — I — 1 2 

;:^i JBiy:4)nN' JB(y:-^)nN' \y - yr 



-C{~eo) I . /^^dy. (2.33) 

iB{y,f)nN' \y-yor 



Thus, by (2.32), (2.33), Lemma 2.6 and a straightforward calculation, we infer that 

/ - ^^^^^dy<C{l + \\PsJ\L^ + \\Pe\Us\'\\L^ + \\u,\\l.) 

Js(g,f)niv' \y-yr 

+C{e) f _ :^^L±^dy + eo||p.u,|Ui+e2||Vp,||2 (2.34) 

JBiyo,f)nN' \y-yr ^ ' 

If we take e sufficiently small, Eq < (5 and e < 6'^, then (2.31) follows from (2.12), (2.13) and 
(2.34). D 

Finally, taking 6o = 6^, we use Lemmas 2.3-2.7 and the finite covering theorem to obtain 
Theorem 2.2. This completes the proof of Theorem 2.2. 

2.3 Vanishing limits as £ — )> and ^o ~^ 

According to [15] and Theorem 2.2, if we take the limits as e — t- and eg ~^ respectively, 
then there is at least a weak solution (ps, u^) to the problem 

dw{psus) = 0, (2.35) 

-fiAus - pVdivus + div{psus (g) u^) + VPsips) = p^f + g, (2.36) 

us\dn = 0, (2.37) 

satisfying (7 = max {7, 4}): 

pseL^^^in), useH^in), I psdx = M; (2.38) 

Jq 

diY[b{ps)us] + [b'{ps)p5 - b{ps)] divu5 = in V'{n); (2.39) 

/ [p\\7us\^ + il\dwus\'^]dx < psf -us+g-usdx, (2.40) 

Jn Jn 

I Ps+ps\^s\\ ^ < ^^^ ^ ll^^ll^^ ^ ||p^|u,|2||^, + ||U5||2,0 (2.41) 

for all a G (0, 1), xq G ^2, where the constant C depends on f, g, /i, /i, M, 7, il and a, but 
not on 6 and xq, 6 is the same as in (1.5). 

Now, with the help of (2.41) we are able to derive the uniform- in-(5 estimates for {ps,us), 
which will be used in passing to the limit as (5 — )■ in the next section, to get a weak solution 
of the system (1.1), (1-2). More precisely, denoting 

^=\\^5\^s\ + Psl^Sl IIli, < a < 1, (2-42) 

we have the following uniform estimates: 



Theorem 2.8. For A defined by (2.4-2), it holds for any 1 < r < 2 — I/7 that 

A + WvisWm + II-P^IIl'- + IIp^Iu^I^IIl'- + IIp^u^IIl'- < C, (2.43) 

5 / p^dx ^0 as (5 ^ 0, (2.44) 

Jo. 

where the constant C depends only on ||f||Loo, ||g||Loo, ^, /i, M, 7 and a (hut not on 6). 

Theorem 2.8 can be obtained by arguments very similar to those used in [12], and hence 
we omit the details of the proof. Instead, we only briefly describe the main steps of the proof 
here. 

Lemma 2.9. Let {ps,us) be the solution of the approximate problem (2.35)-(2.37). Then, 

7-1 

where the constant C depends on ||f||j^oo, ||g||L°°; ^, P, M, 7 and Q, but not on 6. 

Lemma 2.10. Let {ps,us) be the solution of (2.35)-(2.37). Then for s e (l,a + 1 - a/7], 
we have 

\\P5\^Ji5\^\\is + \\P5\\is <c(i + ^^^S^), 

where the constant C depends only on W^Wl^^ , ||g||L°°; P, A, M, 7 and il, but not on 6. 

Lemma 2.11. Let A be defined by (2.42), then we have 

A < C||u5||^i(l + llP^ll^i + ||p5|u5|2||^i + llu^ll^i), 

where the constant C depends on ||f||ioo, ||g||L°°, p, P, M, 7 and 0,, but not on S. 

Remark 2.12. The uniform estimates (2.4 1) is the key estimate when proving Lemma 2.11. 

Recalling that Lemma 2.9 is true for any sG(l,a + l — a/7], we write s = 1 + e, where 
e will be chosen small enough later on, and use Lemma 2.9-Lemma 2.11 to infer that 



A < C||u5|||^i [1 + llu^lll^i + \\ps\us\^\\l^ + WPsh-i 
7— 1 / 7—1 7s— 1 1 \ 

< (7^2(a7+7-2) M _|_ ^2(a7+7-2) _|_ ^ 7a+7-2 ' 1+e J (2.45) 

<C(l+ ^2(7c7+7-2)+'^(=) 



Since (2.45) remains valid for any a G (0, 1), if we write a = 1 — a with < a < 1, then 

1 3(7 - 1) 

^ > =^ — — < 1 

' 1 - 2cr 2(7Q + 7-2) 



where e and a can be arbitrary small. So, by our choice of the parameters e and a, the 

3(7-1) 
2(7Q+7-2) 

A<C, 



exponent 2(^n+^-2) +^(^) ^^^ t)e made less than 1, and therefore we conclude by (2.45) that 



which immediately implies that 
llu^ll/fi < C, 



, ( r 2r\ll2/ r\l/2 ^ r<\\ r 2r||l/2|| r||l/2 ^ ^ 

Jn 
S I pidx < C6^^^&) ( / dp^-^^^'-^Ux^ ^^^^ < CJ^rfe^. 



In ^Jn 

Therefore, we obtain Theorem 2.8. 
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3 Proof of Theorem 1.2 

In this section we will take to the limit as (5 — t- for the approximate problem (2.35)-(2.37) 
to obtain a weak solution of (1.1)-(1.3) for any 7 > 1. The main task for completing the 
proof of Theorem 1.2 is to get the strong convergence of ps to p in L^(r2), and this will be 
fulfilled in several steps. In fact, we have already had all the necessary estimates, and we can 
then employ a process very similar to that used in [15, 12] to get the strong convergence of 
Pi in L^(r2). Hence, we omit the details of the proof and only describe the main steps of the 
proof here. 

First, the following preliminary lemma will be used in the proof of Theorem 1.2. 

Lemma 3.1. ([12]) Let 1 < pi,P2,P < oo, p < pi and Q be a bounded domain inW^. Suppose 
that 

fn^^f weakly in L^^{Q), 
Qn^ g strongly in L^'^iQ), 

and 

fng-n O'fe uniformly bounded in U\Vl). 

Then there is a subsequence of fn9n (still denoted by fng-n), such that 

fnOn -^ fg weakly in ^{Vl). 

Remark 3.1 We point out here that 1 < 1 is not needed in Lemma 3.1. 

^ pi P2 — 

By Theorem 2.8, and the compact imbedding H^{Q,) "^-t-^-)- LP{Q), p G [1,6), we have the 
following limits: 



(3.46) 



It is easy to see that Theorem 2.8, together with Lemma 3.1 and (3.46), yields that 

ps^s -^ pu and psus (X" u^ ^ pu (g) u in L^{Q). (3.47) 

Letting (5 — )■ in (2.35) and (2.37), applying (3.46) and (3.47), we find that the weak limit 
(p, u) of {ps,us) satisfies 

div(pu) = in V'{n), 

div(/9u (g) u) + aVp'y — pAu — /iVdivu = pf in V'{Q); 



6pl^0 


in V'{n), 




Us ^ u 


weakly in Hq{Q), 




us -^u 


strongly in L^{^), 


l<p<6 


ps-^ P 


weakly in L^''{n), 




pJ^T' 


weakly in L^{^). 





dw[b{ps)us] + [b'{ps)p5 - b{ps)]dwus = in P'(i7), 

/ [/i|Vu|^ + /i|divu|^]dx < / pf-u + g-udx. 
Jn Jn 

As the next step, we need to prove an identity (i.e. Lemma 3.2 below) for the so-called 
effective viscous flux defined by 

Hs := apj — {p + p)divus ^ H := ap"/ — {p + /i)divu, as 5 — )■ 0. 
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Lemma 3.2. ([15]) For any (j) G C^{^), we have 

lim / 4){x)HsTk{ps)dx = / 4>{x)HTk{p)dx. 
where Tk{t) is the cut-off function: 

n{t) := A:r(^), T{t) := | *' * ^ ^^' G C^{R), concave. 

Remark 3.2. By using the density argument, one can actually take (/>(x) = 1. 
Lemma 3.3. (Control of the oscillation of the density) We have 

lini||Tfc(/?5) - Tkip)\\Li+^ < C, 
5— s-O 

where the constant C is independent of k. 

Lemma 3.4. (Renormalized continuity equation) The weak limit {p, u) is a renormalized 
solution of (1.1), i.e., {p,u) satisfies 

div[6(p)u] + [b'{p)p - 6(p)]divu = in V'{Vt) 

for any b £ C^(R), h'{z) = for sufficiently large z. 

Now, introducing a family of functions 

{zlogz, < z < k 

GCi(M+)nCO[0,oo), 
zlogk + z Jj^ ^^ds, z>k 

and making use of Lemmas 3.5 and 3.6, we argue, in the same manner as in [1], to conclude 

that 

lim\\ps- pWli =0, 
5— s-O 

which, by (3.46) and the interpolation theory, implies that 

pS ^ p strongly in L^{Q), V 1 < p < 'jr. 

Consequently, we have 

p^ = p^ , a.e. 

Thus, the proof of Theorem 1.2 is complete. 
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